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ABSTRACT 

The intrinsic luminosity functions of extremely fast jets found in many active galaxies and 
gamma-ray bursts are difficult to measure since their apparent luminosities are strongly affected 
by relativistic beaming. Past studies have only provided analytical predictions for the beamed 
characteristics of populations in which all jets have the same Lorentz factor. However, the jets 
found in active galaxies are known to span a large range of speeds. Here we derive analyti- 
cal expressions for the expected Doppler factor distributions and apparent (beamed) luminosity 
functions of randomly-oriented, two-sided jet populations which have bulk Lorentz factors dis- 
tributed according to a simple power-law over a range [Gl, G2]. We find that if a jet population 
has a uniform or inverted Lorentz factor distribution, its Doppler factor distribution and beamed 
luminosity function will be very similar to that of a jet population with Lorentz factors all equal 
to G2. In particular, the beamed and intrinsic luminosity functions will have the same slope at 
high luminosities. In the case of a steeply declining Lorentz factor distribution, the slopes at 
high luminosities will also be identical, but only up to an apparent luminosity equivalent to the 
upper cutoff of the intrinsic luminosity function. Also, at very high apparent luminosities, the 
slope will be proportional to the power law slope of the Lorentz factor distribution. At very low 
luminosities, the form of the apparent luminosity function is very sensitive to the lower cutoff 
and steepness of the Lorentz factor distribution. We discuss how it should be possible to re- 
cover useful information about the intrinsic luminosity functions of relativistic jets by studying 
flux-limited samples selected on the basis of beamed jet emission. 

Subject headings: galaxies : active — galaxies : jets — quasars : general — gamma rays: bursts 
— relativity 



1. INTRODUCTION 

The observed fluxes of highly relativistic jets such as those found in active galactic nuclei and gamma- 
ray bursts are subject to strong relativistic aberration and Doppler beaming effects. This poses a serious 
challenge for determining the intrinsic luminosity functions (LFs) of these objects. Previous studies (e.g., 
Jackson k Wall 1999; Urry & Shafcr 1984; Urry & Padovani 1991; Urry et al. 1991; Urry & Padovani 1995) 
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have discussed in detail how luminosity functions can be apparently altered by Doppler beaming, but present 
analytical expressions only for cases where every jet in the population has the same intrinsic jet speed. In 
their analysis of beamed luminosity functions, Urry et al. (1991) found evidence that BL Lac objects could 
be the beamed counterparts of FR-I galaxies, provided that they have a wide distribution of jet Lorentz 
factors. Large proper motion surveys (Vcrmculcn ct al. 2003; KcUermann et al. 2003) have also revealed 
that the overall AGN population contains a wide range of apparent jet speeds, which is inconsistent with a 
single- valued intrinsic speed distribution (Vermeulen & Cohen 1994; Lister & Marscher 1997). 

Here we extend the earlier work of Urry & Padovani (1991) by deriving full analytical expressions 
for the Doppler factor distributions and beamed luminosity functions of randomly-oriented, two-sided jet 
populations that contain a range of jet speed v. Wc use an improved method that involves first deriving 
the predicted Doppler factor distribution, and then integrating over an appropriate range of Doppler factor 
to obtain the apparent luminosity function. We show how this method results in a much wider number of 
scenarios that can be treated fully analytically. 

For simplicity, we consider only straight jets, where the emission is beamed proportionally according 
to the Lorentz factor of the jet flow (defined as 7 = (1 — (u/c)^)"-'^/^, where c is the speed of light). Wc 
restrict our analysis to cases where the jet Lorentz factors are distributed according to a power law, since 
this distribution provides the best fit to the apparent speed distribution of AGN jets (Lister & Marscher 
1997). Other forms, such as Gaussian and single- valued functions, provide very poor fits to the observed 
speeds. 

In § 2 wc outline our model assumptions and derive the predicted Doppler factor distributions for 
Lorentz factor distributions of the form ex 7". We use these results in § 3 to examine how a simple 
intrinsic power-law luminosity function 4>{£) oc is transformed by Doppler boosting into its observed 
form ^{L). In § 4 we discuss how our results can be applied to recover useful information about the parent 
population of samples dominated by relativistic jet emission, such as gamma-ray bursts, gamma-ray loud 
AGN, and core- dominated, radio-loud AGN. We summarize our findings in § 5. 



2. DOPPLER BEAMING FACTOR DISTRIBUTIONS 

For the purposes of obtaining analytical expressions, we adopt a simple model in which each source 

in the population contains two identical, straight, oppositely-directed jets with bulk Lorentz factor 7. Wc 
assume that the jet axes arc randomly oriented, so that if 6 represents the angle between the jet axis and 
the line of sight, then Pg = sin0, and Q° <6 < 90°, where Pg represents the probability density function of 
the viewing angles. The Lorentz factors have a density function P^ that is positive and non-zero over the 
region 71 < 7 < 72. The kinematic Doppler factor is defined as 

r / 1 -1 

5 = 7-^72 - icos6» , (1) 

and represents the frequency shift between the observer and jet rest frames. The minimum possible Doppler 
factor is Smin — 12^, and occurs when 6 = 90°. The maximum possible Doppler factor in the population is 
^max = 72(1 + /32), where /3 = v/c, and occurs when 9 = 0°. When 72 is large, Smax — 272- 

In Appendix A we derive the expected Doppler factor distributions of randomly-oriented jet populations 
having various Lorentz factor distributions. We plot some individual cases corresponding to single-valued 
and power- law distributions in Figure 1. 
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2.1. Single- Valued Lorentz Factor Distributions 

The two thick lines in Figure 1 represent populations where all jets have the same intrinsic Lorentz 
factor 7o, with the dotted-dashed line corresponding to 70 = 3, and the solid line corresponding to 70 = 30. 
The latter represents the fastest jet speed typically scon in large AGN proper-motion surveys (KcUermann 
et al., in preparation). For a single-valued Lorentz factor distribution, Pg cx 6'^, which corresponds to a 
straight line with slope of —2 in our log-log plot (Fig. 1). 

2.2. Power-Law Lorentz Factor Distributions 

The thin curves in Figure 1 correspond to various jet populations with Lorentz factors distributed 
between 1 < 7 < 30. The thin dotted-dashed line shows the expected Doppler factor distribution for a jet 
population with an inverted 7 distribution (i.e., oc 7). The thin dashed curve represents a uniform 7 
distribution {Pj = constant). These two curves have an overall slope that is roughly the same as in the 
single- valued case (slope = —2) over a wide range of (5, but they both taper off sharply near Smin and Smax- 
Jets with these extreme Doppler factors all have Lorentz factors ~ 72, and viewing angles close to 90° and 
0°, respectively. They are therefore very rare in the population. In the case of models where P^ cx 7" and 
a < 0, the tapering off at the cxtrcma is a result of having relatively fewer jets with 7 = 72, as compared to 
a single-valued population where all the jets have the same jet speed. 

The remaining two curves in Figure 1 represent the cases P^ oc 7^^ (thin solid line) and P^ cx 7^^ (thin 
dotted line). In populations such as these which have few high-speed jets, the distributions taper off much 
more sharply for Doppler factors above 7f ^. 

The slopes of the Pg functions have a discontinuity at 5 = 7f ^, which arises because, if 72 > 71(1 -|- /3i), 
potentially any jet in the population can have a Doppler factor in the range 7^^ < S < 71(1 -|- /3i). This 
leads to a slight over-density of sources in this interval. For a jet to have a Doppler factor that lies outside 
this range, it must have a speed greater than (but not equal to) /3i. 

2.3. Discussion 

We have shown that the power-law slope of the Lorentz factor distribution can have a strong influence on 

the distribution of Doppler factors in a randomly-oriented jet population. In Figure 2 we show the effects of 
increasing the lower limit of the Lorentz factor distribution on a population with 71 < 7 < 72 and P-y oc 7~^. 
As 71 is increased, the region 7^"^ < 5 < + /3i) grows wider, and the cusp at 5 = 7f ^ is replaced by 
a broad region with slope equal to —2. The same behavior also occurs for other values of a. Figure 2 also 
confirms that the Doppler factor distribution approaches that of the single-valued Lorentz factor case as 

71 72- 

It is important to note that in a sample of purely randomly-oriented two-sided jets, exactly half will have 
approaching jet viewing angles greater than 60° . This means that regardless of the form of the jet Lorentz 
factor distribution, there will always be a preponderance of sources with low Doppler factors, such that the 
mean value of 5 will always be less than unity. Somewhat paradoxically, if a jet population is measured to 
have a high mean Doppler factor, this actually implies that it must be dominated by slow jets. We caution, 
however, that this reasoning only applies to genuinely orientation-unbiased samples. Flux-limited samples 
selected on the basis of relativistically-beamed jet emission, on the other hand, will be strongly biased toward 
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Fig. 1. — Predicted Doppler boosting factor distribution functions for randomly-oriented populations of 
two-sided jets. The thick dot-dashed and dotted lines correspond to a jet population having a single bulk 
Lorentz factor of 70 = 3 and 70 — 30, respectively. The thin lines correspond to jet populations whose bulk 
Lorentz factors are distributed over the range 1 < 7 < 30 with a density function cx 7°. The models 
plotted correspond to a = 1 (thin dot-dashed line), a = (thin dashed line), a = — 1 (thin sohd hue), and 
a = — 2 (thin dotted line). 
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small viewing angles, and as a result will always be dominated by intrinsically fast jets with high Doppler 
factors (i.e., blazars; see Vermeulen & Cohen 1994). 



3. BEAMED LUMINOSITY FUNCTIONS 

We now use our results on the predicted Doppler factor distributions to derive the corresponding appar- 
ent luminosity functions of rclativistically beamed jet populations. For the purposes of obtaining analytical 
expressions, we will ignore any luminosity contribution from the jet in each source that is pointing away from 
the observer. In § 3.1 we will show using numerical simulations that for a steep intrinsic LF, the receding 
(counter-jet) luminosity only affects the shape of the faint end of the beamed LF, and has no effect at high 
observed luminosities. 

Due to relativistic aberration, the luminosity of an isotropically-emitting, moving object according to a 
stationary observer is 

L = SPC, (2) 

where C is the source rest-frame luminosity, p = 2 — a for continuous jet emission, and a is the spectral 
index (S^ oc j/"). In the case of a discrete emitting region, we have p = 3 — a due to relativistic time dilation 
of its finite synchrotron-emitting lifetime (Lind & Blandford 1985). For fixed p, the observed luminosity 
function of a population of relativistic jets will depend solely on its Doppler factor distribution and its 
intrinsic (unbeamed) luminosity function. In the figures that follow, we will only plot models with p = 2, 
which is appropriate for the continuous jet emission associated with the flat-spectrum cores of blazars. 

We begin by assuming that the jet population has an intrinsic power-law luminosity function of the 
form: 

fk(£) = / kiJO.~^ for £i < £ < £2 

' 1 , elsewhere, ^ ' 

where the normalization constant is 

-j^ ^ 

^1 = ■ (4) 

The transformation law for probability density functions yields the joint probability function 

hiL,S)=PsiS)0{C)'§^- (5) 
Using equations (2) and (5) and the intrinsic luminosity function, we obtain 

h{L,d) = kiL-^Ps{5) S-P^^-^K (6) 

Integrating over 6 gives 

$(L) = kiL-^ J Ps{d) S^'P+^d6, (7) 
where (using the notation of Urry & Shafer 1984) 

C = B- (l/p) - 1 . (8) 

In their derivation of the beamed LF for a jet population having a unique Lorentz factor, Urry & Shafer 
(1984) identified several important values of L. Here we keep their definitions of L3 and L4, and extend their 
notation as follows: 




Doppler Factor (6) 

Fig. 2. — Predicted distribution functions of Doppler factor for jet populations with bulk Lorentz factors 
distributed according to (x 7^^, and ranging up to 72 — 30. The solid curve shows the predicted 
distribution for a lower limits of 71 = 1, while the dotted and dot-dashed curves represent models with 
71 = 3 and 71 — 10, respectively. 
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Wc note that a similar extension of notation was made by Urry & Padovani (1991), however, they do 
not include the parameter L5, which is needed for the analytical solutions presented in Appendix B. 

When 72 is large, the beamed LF spans a wide range: 



Thus, even a jet population with a single intrinsic luminosity £ = £1 = £2 will end up with a broad beamed 
luminosity function. This scenario has been invoked to explain the wide range of observed luminosities in 
gamma-ray bursts, which are thought to have Lorentz factors on the order of 100 (e.g., Kumar & Piran 



To determine an analytical expression for the beamed LF, we substitute the relevant expression for 
Ps (see eq. A9) into equation (7) and integrate over appropriate limits. We find that analytical solutions 
are generally only possible when the product Cp is an integer. We derive expressions for several cases in 
Appendix B. For simplicity, we have considered only cases where B > {p + l)/p and L3 > L4. This 
corresponds to the condition {C^/jCiY^p > 72(1 + (32), which will always be satisfied provided £2 is made 
arbitrarily large. The latter can be justified if, as in our case, the intrinsic LF slope is fairly steep (i.e., 
B > 1.5 for p = 2). With these choices of parameters, the L values in equation (9) are arranged in order of 
increasing apparent luminosity. 

In Figure 3 we plot log [L ^{L)] versus log[L/£i] for various power-law Lorentz factor distributions 
distributed over the range 1 < 7 < 30. The thick solid curve is the original (unbeamed) LF, which in this 
example has a power-law slope equal to —2 and log [£2/£i] = 7. The remaining curves represent the beamed 
LFs for various vahies of power-law index a, where oc 7". 

Urry & Sliafer (1984) found that for a single- valued Lorentz factor distribution, the beamed luminosity 
function will have the same slope as the unbeamed one for L > L4, except in the immediate vicinity of Lmax- 
In their numerical analysis, Urry & Padovani (1991) claimed that this also holds true for a distribution of 
Lorentz factors. We confirm this in the case of flat or inverted Lorentz factor distributions (i.e., o = and 
a = 1). Indeed, the beamed LFs of a single- valued model with 7 = 7o and a uniformly distributed model 
(P-y oc 7 ; 1 < 7 < 7o) are nearly identical for L > L4 (see Fig. 4). 

For steep power- law 7 distributions with a < — 1, the beamed and unbeamed LFs also have the same 
slope, however, this is true only for the region L4 < L < Lg. Above Lg, the slope steepens by an amount 
proportional to the power-law 7 index a. We note that in the case of a jet population where C2/ Ci is large 
and the intrinsic LF is steep, objects with L > Lg will be extremely rare. In the models shown in Figure 3, 

only one jet in ^ 10^ will have L > Lg. 

At luminosities fainter than L4, Urry & Padovani (1991) claimed that the slope will always be equal 
to —{p + l)/p, regardless of the Lorentz factor distribution. However, we find this is true only in cases 
where there are very few low-Lorentz factor jets in the population (i.e., when 71 >> 1 or a > 0). In cases 
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where 71 ~ 1 and the Lorentz factor distribution is steep, the LF slope at low luminosities is significantly 
flatter than — (p+ l)/p, and can even be highly inverted (see Fig 3). We caution, however, that at these 
low-luminosities many jets will lie in the plane of the sky, and that our analytical expressions ignore the 
receding jet. Adding the counterjet's contribution to the total emission will change the predicted luminosity 
functions slightly. We discuss this issue in more detail in the next section. 
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Fig. 3. — Predicted apparent luminosity functions (LFs) for jet populations with Lorentz factors distributed 
according to a simple power law. The intrinsic LF (thick solid curve) has a power-law slope of —2, and 
ranges from Ci to £2, with C2/ Ci — 10^. The remaining curves represent the beamed (apparent) LFs for jet 
populations with Lorenz factors distributed according to cx 7" over the range 1 < 7 < 30. These models 
plotted correspond to a = 1 (thin sohd curve), a = (dotted curve), a = — 1 (dot-dashed curve), and a — —2 
(dashed curve). The specific locations of L3, L4, and Lg are marked with arrows (see eq. 9). 
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Fig. 4. — Comparison of beamed luminosity functions for jet populations with single-valued and power-law 
Lorentz factor distributions. The thick solid line represents the intrinsic LF, which has a power-law slope of 
—2, and ranges over £i to £2, with C2/ Ci = 10^. The other curves represent beamed LFs for a jet population 
with a single Lorentz factor 70 = 30 (thin solid curve), one with Lorentz factors uniformly distributed over 
the range 1 < 7 < 30 (dotted curve). 
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Fig. 5. — Comparison of beamed luminosity functions for jet populations with Lorentz factors uniformly 
distributed over different ranges (i.e., P.y oc constant, where 71 < 7 < 72)- The thick solid line represents 
the intrinsic LF, which has a power-law slope of —2, and ranges over Ci to £2, with C2/C1 = 10^. The 
other curves represent the beamed LFs of jet populations with 1 < 7 < 10 (thin solid curve), 1 < 7 < 30 
(dot-dashed curve), 20 < 7 < 30 (dotted curve). 
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We have also investigated the effects of various ranges of Lorentz factor on the apparent luminosity 
functions of jet populations. In Figure 5 we plot three different models with Lorentz factors uniformly 
distributed over different ranges 71 < 7 < 72- It is apparent that for L > L4, the bright end of the beamed 
LF is largely insensitive to the values of 71 and 72. Increasing 72 can potentially increase the number of 
ultra-high luminosity objects, but again these are likely to be extremely rare for steep intrinsic LFs. 

At faint luminosities, the beamed LF is highly sensitive to the value of 72, since this dictates the 
minimum Dopplcr factor, and hence L^nin- Because a very fast jet population will have a low mean Doppler 
factor (see § 2.3), its sources will, on average, appear fainter than those of a population dominated by slow 
speeds. Again, this reasoning applies only to genuinely orientation-unbiased samples. 

3.1. Two-sided jets 

The beamed luminosity functions presented in the previous section are not strictly correct at low lumi- 
nosities, since they ignore any luminosity contribution from the receding jet in each source. When considering 
only highly beamed objects such as blazars or gamma-ray bursts, this is a reasonable approximation, since 
the receding jet will be extremely faint (and usually undetectable). However, in samples containing mildly 
relativistic jets, or jets that lie close to the plane of the sky (e.g., lobe-dominated radio galaxies), both the 
approaching and receding jets can contribute roughly equal amounts of flux. 

In Figure 6 we compare two models with B = 2, a = —1, 1 < 7 < 30, and £2/^1 = lO*". The solid hue 
shows the beamed luminosity function if only the luminosity of the approaching jet is considered. The dashed 
line represents a numerically-determined beamed LF in which the luminosities of both the approaching and 
receding jets are taken into account. For luminosities below Li, the curves differ considerably, since these 
jets lie close to the plane of the sky. At higher luminosities, the predicted LFs are virtually identical, since 
the apparent approaching/receding jet luminosity ratio becomes very large for small viewing angles and 
high Lorentz factors. We find that this result holds true for all ranges of B and a discussed in this paper. 
We conclude therefore, that the emission from the receding jet can be safely ignored when studying the 
luminosity functions of highly beamed objects. 

4. OBSERVATIONAL IMPLICATIONS 

The analytical results derived in this paper are applicable to any astronomical sample that is selected 
solely on the basis of relativistically-beamed jet emission. Our results are particularly relevant to compact 
radio- loud AGN, whose high-frequency radio emission originates primarily from milliarcsecond-scale rela- 
tivistic jets. With the advent of dedicated VLBI arrays, it is now possible to conduct large multi-epoch 
surveys of these objects (e.g., the Caltech-Jodrell Surveys; Taylor et al. 1996 and the VLBA 2 cm Survey; 
Kellermann et al. 1998). However, because of the high thresholds required for fringe detection, these large 
VLBI studies have so far been confined to the brightest, most compact objects whose radiation is highly 
beamed toward us. These types of samples are not ideal for studying the general properties of the AGN jets, 
since one miist first disentangle the combined effects of beaming, jet axis orientation bias, and Malmquist 
bias before one can construct a luminosity function or intrinsic speed distribution. Monte Carlo simulations 
have played an important role in this respect. Lister & Marscher (1997) showed that based on their observed 
superluminal speed distribution, AGN jets must have a wide range of Lorentz factors distributed according 
to a relatively steep power law {P^ oc j"', where a ~ —1.5). They also demonstrated that a reasonably 
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Fig. 6. — Beamed luminosity functions for two relativistic two-sided jet populations with identical Lorentz 
factor distributions and intrinsic luminosity functions. The dashed line is calculated numerically, and shows 
the beamed LF if the luminosity from both the approaching and receding jets in each source are considered. 
The solid line shows the expected LF if the receding jet is ignored. The two curves differ only at low 
apparent luminosities, where the Doppler factors are small and the approaching and receding jets have 
similar luminosities. 
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large, flux-limited sample of bright jets should always contain some objects with apparent speed equal to 
the maximum jet Lorentz factor in the parent population. 

Since the observed LF of radio-loud AGN is fairly steep (B ~ 1.9±0.1; Willott ct al. 1998), our findings 
suggest that the observed LF of a complete flux-limited sample selected on the basis of beamed emission will 
closely match the intrinsic LF for luminosities in the range 4£i/3^ ^^^^ ^ L ^ £2- Here we have assumed a 
boost index p = 2 and 72 equal to the maximum apparent speed Pa.max- The lower luminosity limit £1 can 
be constrained using the nearest sources in a complete VLBI flux-limitcd sample, since these will generally 
meet the flux density criteria on the basis of their cosmological proximity, and not their Doppler boosting 
factors. The upper luminosity limit C2 can be constrained using the observed space density of radio-loud 
AGN. 

We are currently gathering VLBI proper motion and flux density data at 15 GHz on a sample of the 
brightest ^ 130 jets in the northern sky. This sample, entitled MOJAVE (Monitoring Of Jets in AGN with 
VLBA Experiments; Zensus et al. 2003), is complete on the basis of compact (VLBI) flux density, and is 
therefore well-suited for exploring the intrinsic luminosity function of AGN jets. An analysis of the MOJAVE 
data is currently underway, and the results will be discussed in an upcoming paper. 

5. SUMMARY 

Astronomical samples of objects selected on the basis of emission from relativistic jets, such as gamma- 
ray bursters, gamma-ray loud AGN, and compact radio-loud AGN, will by highly affected by relativistic 
beaming effects. In particular, their observed luminosities are predicted to be c^nhanccd by a factor 5^ , 
where 5 is the Doppler factor. We have derived expressions for the predicted Doppler factor distributions 
of two-sided, randomly-oriented, straight jet populations whose bulk Lorentz factors are distributed in the 
range [71,72] according to a simple power-law oc 7". We summarize our flndings as follows: 

• The minimum and maximum jet Doppler factor in the population will depend only on the maximum 

Lorentz factor 72. Speciflcally, 5min = 12^, and 5max = 72(1 + /32)- 

• An orientation-unbiased sample of jets will always have a mean Doppler factor less than unity, since 
the majority of jets will lie close to the plane of the sky. Accordingly, a jet population weighted toward 
very fast speeds will have a lower mean Doppler factor, and therefore appear fainter than an identical 
jet population with slower intrinsic jet speeds. 

• A randomly-oriented jet population with a flat (a = 0) or inverted (a = 1) Lorentz factor distribution 
will have a Doppler factor distribution very similar to a population where all jets have 7 = 72- Over 
a broad range of 5, the distribution function will vary as (5^^, except near 5min and dmax) where 
much fewer sources are expected compared to the single- valued case. Sources with the highest possible 
Doppler factors near 5max will always be rare, even when all the jets in the population are highly 
relativistic (i.e., 71 >> 1). This is a consequence of the extremely small angles to the line of sight that 
are necessary for high Doppler factors. 

We have also derived analytical expressions for the apparent luminosity functions (LFs) of relativistic 
jet populations having a power-law distribution of Lorentz factors and a steep intrinsic power-law luminosity 
function. We find that: 
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• The slope of the apparent LF at low luminosities is not always equal — (p + as claimed by Urry 
& Padovani (1991). The apparent LF in this region is highly sensitive to the speed distribution of the 
parent population. In particular, if the parent population is dominated by slow jets, the LF slope will 
be very flat or inverted at low luminosities. 

• Regardless of the values of 71 and 72, when the Lorentz factor distribution of the jet population is 
steep (i.e., a < —1, where oc 7"), the intrinsic and beamed LFs will have the same slope over the 
range LiS^^^^ < L < C2, where L is the observed luminosity and £1 and £2 are the lower and upper 

cutoffs on the intrinsic LF, respectively. 

• Since the intrinsic LF of radio-loud AGN jets is known to be steep, it should be possible to investigate 
its properties using flux-limited VLBI samples of highly beamed AGN jets. Such samples will have 
to be selected at relatively high radio frequencies (e.g. above ^ 10 GHz) in order to eliminate any 
contribution from un-beamed, steep-spectrum emission. An analysis of such a sample selected at 15 
GHz (MOJAVE) is underway, and will be presented in an upcoming paper. 
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A. DOPPLER FACTOR DISTRIBUTIONS OF JET POPULATIONS WITH A POWER 

LAW DISTRIBUTION OF LORENTZ FACTORS 

We wish to determine the probability density function Ps{S), which describes the expected distribution of 
Doppler factors for a randomly-oriented, two-sided jet population. The theory of probability transformation 
for several variables implies 

9{d,'y) = P^{'y)Pe{0) '^^ 



d6 



Now 



5V7^ - 1 sin6». 



Since the viewing angles are distributed according to 

Pg{0) = sine, 

we have 



We integrate over 7 to obtain the probability density function for 5: 



J f 



f{5) - 1 



(Al) 
(A2) 

(A3) 

(A4) 

(A5) 



where the limits on this integral are illustrated graphically in Figure 7. The lower limit depends on the value 
of S as follows: 

for 72-1 <S < 7f 1 

71 for 7r^ < < 71(1 + /5i) 
iiip for + f3i) < 6 < 6max 

72 elsewhere, 



m = 



(A6) 
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where we have used the fact that -\/7^ — 1 = 7/3 and 5max = 72(1 + 02)- 

Analytical solutions of equation (A5) arc generally possible if oc 7", where a is an integer. Below we 
derive expressions for Ps{S) for power-law slopes ranging from a = -|-ltoa = — 2. 



A.l. Case 1: Single- valued Lorentz factor distribution 

When 71 = 72 = 7o, then P-yij) can be considered as a delta- function centered at 7 = 70. Equation (A5) 
yields the well-known result (e.g., Urry & Shafer 1984): 



Ps{S) 



S'^ lo' Po' for lo' < ^ < 7o(l + /3o) 
, elsewhere. 



(A7) 



A. 2. Case 2: Uniform Lorentz factor distribution 

In the case where the Lorentz factors in the jet population are uniformly distributed within the range 
71 < 7 < 72, this is equivalent to setting a = 0, i.e.: 



^7(7) 



_/(72-7i) ^ for 71 < 7 < 



72 



elsewhere, 
we substitute this function into equation (A5) and perform the integration to obtain 

for 72^1 <S< 7f 1 



1 



In 



^2 - 71) 



X < 



In 
In 
vO, 



for 7r' <'5<7i(l+/3i) 



7i(l+/3i 
'-'max 

/5\ for ii{l + (3i)<5<5 

r? 

elsewhere. 



(A8) 



(A9) 



For large 72 and 71 ~ 1, equation (A9) reduces to: 

2(572 



1 ^ J In 
•5^72 I In 



for 72'<^<1 
272 /(5| for l<6< 272. 



(AlO) 



A. 3. Case 3: Power Law Lorentz Factor Distributions 
A. 3.1. Power law distribution of the form ^7(7) oc 7 
For the case a = +1 we have 



lo 



for 71 < 7 < 72 
elsewhere. 



(All) 




Doppler factor (5) 



Fig. 7. — Plot of Lorentz factor versus Doppler factor illustrating integration limits for determining (see 
equation A5). For a particular value of 5, the integration limits on 7 are defined by the polygonal region 
(i.e., the lower integration limit is equal to in region A, 71 in region B, and (1 + d'^)/2S in region C.) 
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B. BEAMED LUMINOSITY FUNCTIONS OF JET POPULATIONS WITH A 
POWER-LAW DISTRIBUTION OF LORENTZ FACTORS 



For 7i ^ 1 , the form of the beamed luminosity function is: 



F[(L/C^Y''']-F['t:^\ 

HliL/Ciy'"] + F[7-i] - _F[7-i] + G[7i + 7i/3i] - G^r'] " Hbi + 7i/3i], 
= fciZ,-^ X { F[7-i] - F[7-i] + G[7i + "nPi] - Gh^^] + -H"[<5™ax] - ff[7i + 7i/9i], 

-F[(L/£2)'^''] + F[7-i] + G[7i + 7i/3il - Gh"^] + ff[i5„„x] - //[71 + 7i/3i 
-G[(L/£2)'-'''] + G[7i + 7i/3i] + if[5ma«] - H[-y^ + 7i/3i], 
-H[(L/£2)'''''l + 
, 



Lrn.i-a < L < 

L5 < L < Le 
Le, < L < Li 
Li < L < Ls 
L3 < L < L7 
Lt < L < Lg 
Ls< L< L„„^ 

L > Lmax- 



(Bl) 



If 7i = 1, then = Lq = Ci, and = Lg = Li, and equation (Bl) reduces to: 



X < 



0, 

J^[(L/A)Vp]-F[72-^], 
ff[(L/A)Vf]-F[l]+i^[l]-F[7-i] 

i^[l]-i^[72-']+i?[5™a.]~iy[l], 

[(L//:2)i/p] + F[i] + i/[(5„,,] - £r[i], 

-if[(L/£2)'/P]+i?[(5maJ, 

0, 



Lmin ^ L < Ci 

Ci<L<L4 
L4 < L < L3 

Ls,<L < £2 

£2 ^ L < Ljnax 



(B2) 



The functions F(x), G(x), and -ff(x) are tabulated in the following sections for various integer values 
of a and Cp. They are grouped by particular values of a, where -P7(7) oc 7" in the jet population, and are 
further subdivided by particular values of the product Cp (see eq. 8). 



B.l. Power-Law Lorentz Factor Distribution of the form ^7(7) oc 7 

When a = +1 we have: 



F{x) 



X < 



arcsm x 



In |a;| 

72/322; — Vl — — In 
1 



\l2P2x'^ — \ arcsinx — ^Vl — ■ 

[\l2f32X^ + \{l-X^f/\ 



if 


Cp = 





if 


Cp = 


1 


if 


Cp = 


2 


if 


Cp = 


3, 



(B3) 



G{x) 



2(72/32-7i/3i) flnl'^l. if 



(7l 



■7f) 



^ if 





1,2,3, 



(B4) 



H{x) = 



72 -7i 



72/32 In |a;|- f - ^, 
72P2X - ^ - i In Ix 

Gp 



2(Gjj+l) 2(Cp-l)' 



if Cp = 

if Cp = 1 
if Cp = 2,3. 



(B5) 
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B.2. Uniform Lorentz factor distribution of the form Pj{'y) = constant 
In the case of a uniform 7 distribution (a = 0), we have 



: In 



F{x) = 



1 



72 - 71 



In 



arcsmx, 



— — i arcsinx + 3- In 



if Cp=l 
if Cp = 2 
if Cp = 3, 



G(x) = 



c<^Pln 



7i(l+/3i) 



(72 - 7i)C'P 



if Cp=l,2,3, 

if Cp=l,2,3. 



(72-7i)(Cp)2 

There are no analytical expressions possible when a = and Cp = 



(B6) 

(B7) 

(B8) 



B.3. Power-Law Lorentz Factor Distribution of the form ^7(7) oc 7 ^ 

When a = —1, we have 

arctan (^^) - arctan (^) 



F{x) 



In 1 72/71 



arctan I 



arctan 



(vT^) 



if Cp=l 
f VI -x'^ - I arcsin [x] , if Cp = 2 (B9) 



G{x) = 



Cp In 1 72/71 



arctan 



/ 1 



V7iA 



arctan 



1 



72/^2 



if Cp = 3, 
if Cp=l,2,3, (BIO) 



H{x) 



1 



In I72/71 



arctan 
arctan 
arctan 



\x^-l ) 



2x 

7^ 



arctan 

— arctan 

— arctan 



2/32 



V72/ 

V 72/32 J 



+ ln|a;2 + 1|, 
+ X — arctan a;, 



+ ^-|ln|a;2 + l| 



There are no analytical expressions possible when a = — 1 and Cp = 0. 



if Cp=l 
if Cp = 2 
if Cp = 3. 



(Bll) 



B.4. Power-Law Lorentz Factor Distribution of the form ^7(7) cx 7 
For the case a = —2, we have 



F{x) 



1 



7i 



72 



I arcsin a; - |(2a;^ - - x'^, 



(/32- l)lna;- Vl-a;2 + ln|l 
02X — f a/I — x"^ — ^ arcsin X, 

2 
3 



if Cp = 

if Cp = 1 

if Cp = 2 

if Cp = 3, 



(B12) 
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G{x) 



(/32 - 



X 



{ 



ln|a;|, if Cp = 

if Cp= 1,2,3, 



(B13) 



H{x) = 



7^-72 



1 



-1 



X 



(/32 + l)lna;-ln|l + a;2|, if Cp = 

(/?2 ^ 1)2; + 2 arctana;, if Cp = 1 

i(/32-l)x2 + ln|l+a;2|, if Cp = 2 

i (/32 - + 2a; - 2 arctan a;, if Cp = 3. 



(B14) 



REFERENCES 



Jackson, C. A. & Wall, J. V. 1999, MNRAS, 304, 160 

Kellermann, K. I., Vermeulen, R. C, Zensus, J. A., & Cohen, M. H. 1998, AJ, 115, 1295 

Kellermann, K. I. et al, 2003, in preparation 

Kumar, R & Piran, T. 2000, ApJ, 535, 152 

Lind, K. R. & Blandford, R. D. 1985, ApJ, 295, 358 

Lister, M. L. & Marsclicr, A. P. 1997, ApJ, 476, 572 

Taylor, G. B., Vermeulen, R. C, Readhead, A. C. S., Pearson, T. J., Henstock, D. R., & Wilkinson, P. N. 
1996, ApJS, 107, 37 

Urry, C. M. & Padovani, P. 1995, PASP, 107, 803 

Urry, C. M. & Padovani, P. 1991, ApJ, 371, 60 

Urry, C. M., Padovani, P., & Stickel, M. 1991, ApJ, 382, 501 

Urry, C. M. & Shafer, R. A. 1984, ApJ, 280, 569 (US84) 

Vermeulen, R. C. & Cohen, M. H. 1994, ApJ, 430, 467 

Vermeulen, R. C, & Britzen, S. 2003, in Radio Astronomy at the Fringe, Proceedings of a Workshop on 
Extragalactic Radio Sources held in Green Bank, WV, Oct 10-12, 2002, eds. J. A. Zensus, E. Ros, & 
M.H.C. Cohen, in press. 

Willott, C. J., Rawlings, S., Blundell, K. M., & Lacy, M. 1998, MNRAS, 300, 625 

Zensus, J. A., et al. 2003, in Radio Astronomy at the Fringe, Proceedings of a Workshop on Extragalactic 
Radio Sources held in Green Bank, WV, Oct 10-12, 2002, eds. J.A. Zensus, E. Ros, & M.H.C. Cohen, 
in press. 



This preprint was prepared with the AAS lATJiX macros v5.0. 



